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SUMMARY 


The report is mainly concerned vith the principal investigator's 
research on some analytical problems in the interactions between the 
mean-shear flows and the acoustic field in the planar and circular jets. 

These problems are basic in understanding the effects of coherent large 
structure on the generation and complications of sound in a sub-sonic jet. 

Three problems have been investigated. The first problem pertains to a 
spatial (vs. temporal) normal-mode analysis in a planar jet. This basic 
problem is formulated as an eigenvalue -value problem. Since it is not of 
Sturm-Lioville type, the wave numbers (eigen-values) behave in a much more 
complex manner. This difficulty dictates the use of various methods of 
approximation. Here we adopted an asymptotic method, known as the WKBJ 
method. Eventhough the method was designed for the high frequency case, 
it is often found useful also for lower frequencies, as a result of analytic 
continuation. In reality the disturbance consists of many different frequencies 
Therfore one must deal with the wave packets. Given the spectral function 
of the source, we showed how to determine the excited wave packet by the 
Fourier-method and the method of stationary phase. In the second problem, 
we consider a slightly divergent, planar jet. Since, from a physical viewpoint 
the exponential-growing wave violates the radiation condition. The unstable 
waves must have a bounded envelope in space. In applied mathematics, this 
defect is a manisfestation of a non-uniform perturbation analysis. To obtain 
a uniform expansion, we recommend the two-variable method. By applying this 
method, we derived the evelope equation for the excited waves, which is 
amenable to numerical solution. Thereby the effects of divergent shear flow 




on the sound amplitude can be assessed. Then we extend our analysis to a 

parallel or non-parellel cylindrical jet. The third problem is concerned 

with the acoustic waves in an axisymmetrical jet. By using the cylindrical 

coordinates, we have shown that the analytical techniques used in treating 

the planar jet, such as the WKBJ method, the method of stationary phase and 

the two-variable method, can be applied to the parallel or non-parallel cylindrical 

jets as well. 

In the appendix, a recent paper entitled "Reconstruction of ‘■he Mutual 
Coherence Functions for a Moving Source", is attached. This paper summarize 
a related work partially supported by this grant, on determining the acoustic 
source structure in a moving jet. The method is based the high dimensional 
radon transform. Also it was shown that, in the case of line source, the 
problem becomes similar to that in the X-ray tomagraphy. 
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Introduction 

Recently many workers in jet-noise research have investigated the role of 
the shear mean- flow in the generation of sound. In supersonic flow, the modified 
Lighthill's equations had been derived by Phillips [l], Lilley [2] and Pao [3], 
to account for the shear- flow term in sound generation. The so-called shear 
noise has also been discussed by Ribner [4], Mani [5] and others. All of 
these studies are qualitative in nature. The experimental results of 
Crow-Champagne [6] has stimulated a great deal of interest in the effects of 
coherent shear structure on the acoustic amplification in the jet. In this 
[ paper we shall be concerned with some basic problems in the acoustic-shear flow 

r 

interaction. 

Perhaps, the simplest model for acoustics-flow interaction problem with 
physical meaning is acoustic wave propagation in a two-dimensional parallel 
flow of an inviscid, incompressible and non-conducting fluid. This problem 
is similar to that of stability of two-dimensional parallel flows. A quick 
reading of the recent book on hydrodynamic stability (Chapter 4) by Drazin 
and Reid [7] would reveal that, even in this simple situation, how little is 
known about the general properties of its solution. This impression has motivated 
us to re-examine the two-dimensional model. A number of authors had addressed to 
the growth of spatial waves in parallel flows. However, due to radiation condition 
at « , a spatially 

* This work was supported by the NASA/Langley Research Center under the 
grant NSG - 1330. 
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growing plane wave in a parallel flow seems physically meaningless. A good 
discussion of this point may be found in [7]. It seems more reasonable to 
account for the spatially growing modes by the divergence of shear flows. 

This is the subject of several studies, e.g. [8] - [13]. 

This paper constitutes a summary of our results concerning some acoustics- 
flow interaction problems. In section 1, the governing equations are presented; 
Sections 2 and 3 are concerned with acoustic waves in a planar jet of constant 
width. Section 4 pertains to a slowly divergent planar jet. Here we shall 
introduce a two-variable expansion scheme to derive an evelope equation for 
the spatial growth (or decay) due to the divergent flow. Finally, in Section 5, 
we briefly discuss how to generalize our results for planar jets to the 
axisymmetrical ones. 


1. Governing Equations for the Mean Flow. 

Consider the effect of mean flow on the sound propagation in an inviscid, 
non-heat-conducting jet flow. Let U,u denote the mean and acoustic velocities; 
p o ,p the mean and acoustic pressures and p Q ,p the mean and perturbed densities, 
respectively. By perturbing the Euler's equation, the continuity and the energy 
equation and assuming that the mean pressure gradient terms are small, the 
following set of acoustic equations may be obtained, 


jjf + yp 0 (V • u) = q , (1.1) 

+ T“ 7p + (u • V) v =0 . (1.2) 

uc P 0 


2 


where y 
and 


c /c is the specific heat ratio, q is an applied acoustic source, 
P v 


_D_ 

Dt 



+ v • 7) 


(1.3) 


stands for the material derivative. 

For convenience we use the notations U = (U,V,W), u * (u,v,w) and 
■x * (x,y,z) or (x,r,0) , depending on the rectangular or cylindrical coordinates 
arc used. The following special cases will be considered 

(A) Two dimensional jet: W = w = z * 0 

(i) Parallel flow: V - 0, U - U(y) , 

(ii) Divergent flow: V « V(x,y), U * U(x,y) . 

(B) Cylindrical jet: x = (x,r,0) 

(i) Parallel flow: V = W - 0, U = U(r) , 

(ii) Divergent flow: W ■ 0 , U = U(x,r), V * V(x,r) . 


In what follows, we shall give a more detailed description of Case A, 
while Case B will be discussed only briefly. 


2. Spatial Modal Analysis in a Planar Jet. 

Let us first consider the most basic problem: the normal mode analysis in a 

two-dimensional jet. In components, the governing equations (1.1) - (1.2) 
with q = 0 read 
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hr + 'vS + ^ ■ 0 • 

°u + ^-|^ + VJJ - 0 , -»<*<«, *< y <b 


(2.1) 


Dv + * 0 , 

p c a y 


which is subject to appropriate boundary conditions. 


where 


D ‘<tt + V &' " • "w ■ 


( 2 . 2 ) 


According to the normal-mode stability analysis [ 7 ] , we seek a solution of 
(2.1) in the form for a < y < b 

(P.u.v) » (p,u,v)(y)e i ^ kx “ tut ^ . (2.3) 


For a given frequency uu , we are interested in the dependence of the wave number 
k on uu . If k(uj) is real, the mode is propagating or proper. Otherwise, 
k is complex and the corresponding mode is non-propagating. If Im k a 0 , the 
mode decays and, hence, stable. However, if Im k < 0 , the mode grows 
exponentially as x -* ® . But this means, for fixed t, the modal amplitude will 
grow along the downstream. It is physically meaningless to label such mode as 
unstable, since it violates the boundary conditions at as . To account for wave 
behavior, instead of simple mode (2.3), we shall later replace it by a wave 
packet by superposing modes of different frequencies. 

Let us substitute (2.3) into (2.1) to get a system of ordinary differential 
equations in y 


ikT|p +yp o (iku + ^) - 0 , 


ikT^u + T|'v + ikp/p =0, a<y<b 


ikTlv + “ j * 0 , 

p 0 dy 


(2.4) 
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where 

1] * U(y) - c, with c - (uu/k) . 

After eliminating u,v in favor of t> , we obtain a second order equation 
T)4" - 2Y|'4>’ - k 2 11<l " 0, a<y <b 

where 

4 " p/p d . 

li - Va , or - (YP 0 /P 0 ) 1/,Z 

and 

( “ dy • 

By a simple change of variable 

* - Tit . (2.9) 

the equation (2.6) may be rewritten as 

L (k) t|r * ( 1 - jj^) t|r - rt|f * 0, a<y<b , 

where 

r - lf 2 (T -271' 2 ) . 

Here we assume 1) ■ U(y) - c j 4 0 in (a,b) . Otherwise we let c be complex 
and define \|r as the principal-value as Imc-»0. Consider the homogeneous 
BVP (boundary-value problem) associated with the original system (2.1) 

L(k)f - 0 , 

B^ - B 2 f - 0 , 

where (2.13) stands for appropriate B.C.'s (boundary conditions). Let the 


i.AlailWM 


eigenfunctions f with eigenvalues k R , form a complete set in l/(a,b) , 
the space of square-integrable functions over (a,b) . Suppose we multiply 
(2.12) by i|t , the complex conjugate of i|t , and integrate from a to b 
to obtain 

k 2 J* b [l- (^) 2 ]l*l 2 dy + 2k(%f b U j i|( 1 2 dy + / b {l*’| 2 +tr(k>(“) 2 ] l*| 2 }dy - 0 . (2.14) 

a « */■ a a “ 


The difference between (2.14) and its conjugate yields 

(k 2 -k* 2 )f ll-(^^-) 2 ] | i|;(y) | 2 dy +2(k-k*)r (jfj)U(y) | i|r(y) 1 2 dy - 0 . (2.15) 

w a « a a 

The above show that the eigenvalues need not be real, even though L(k) is self- 
adjoint. Note that (2.14) may be written as 


b k + 2b. k + b » 0 , 
0 i z 


(2.16) 


where , b, , b„ are the corresponding coefficients in (2.14). Thus, if r i& 
o 1 z 

independent of k , corresponding to each eigenfunction i|r n , there are two eigenvalues 

k which are roots of (2.16) . Let 
n 


A(*) « (b‘ - b o .b 2 ) 


(2.17) 


be the discriminant for (2.16). Thus, if A(ijt) * 0 , the roots are real. The 
condition for n>th mode being proper at a fixed uu is that 


M^ n ) * 0 , 


otherwise k Q are complex so that this mode either grows or decays in x 


As a simple example, consider the BVP 
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L(k)| ■ 0, -b*a<y<b , 
f(-b) - *(b) - 0 , 


( 2 . 


[ 

I 


2 2 

where U is constant. Then L(k)i|> ■ V' • k (1 - p )| , and the eigenfunctions are 
{Sin ^(y+b)} * corresponding to eigenvalues k R determined by (2.14) 

[1 - (£) 2 lk 2 + 2(-^)u»k + [® 2 - (Jj) 2 ] - 0 (2. 

Of 

for which 


A/, \ ,U\2 

A(i|r n ) - (“) 


[1 -<S> 2 l<2b 


mi) 2 


1 , 2 ,. 


Thus , for a fixed u) , there are a finite number of proper modes k fl with n S N 
where N is the largest integer less or equal to ■ — with m ■ — , 

emd-m 2 )* « 

The rest of the modes are improper and may grow or decay. 

Next, suppose the velocity profile U(y) is symmetric and varies smoothly 
over [-b,b] . As shown by (2.19), k - 0(n) for large n , the behavior of 
higher-order modes may be determined by an asymptotic method. The well-known 
WKBJ method can be applied. Let 

f ~ A(y)e k ®^ as k -» • , (2.2( 

where A, 0 are determined by the d.e. (2.10). The general solution can then be 
written as 

♦ ~ (l-u 2 )” 1/4 {CSink 9 + C 2 CoshkO) (2.2] 

where 

0 - J y [l-U 2 (s)] l/2 ds . (2.22 

a 

The approximate solution is valid only when 

H> ■ Of' l (u - c) + +1 , a <y < b . (2.2: 


The points y at which j*(y ) * +1 are called turning points. In this case 


t<mm* ** I MftU ■ 
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Che aoluClon (2.21) is not valid near y ■ y Q , and a different form of solution 
must be constructed. Introduce an appropriate change of variables [14] 

§ - f(y) , C - g(y) *(y) , (2.24) 


the equation (2.3) yields 


*"§ - k 2 SC ‘ aC , 
dV 


(2.25) 


where 


5 - f 


[| J y (l-H 2 ) 1/2 dy'] 2/3 , 


g - If'(y)l 


1/2 


(2.26) 


« l til . 1 LH12Z . 


2 (f ) 3 


4 (f ) 4 


(f ) 2 


For large k , the equation (2.25) is approximately an Airy equation so that 


C~ CA 1 (?) + dB 1 (5) as |k| -» 


(2.27) 


where A^,B^ are Airy functions of first and second kind. Hence the solution to 
(2.3) is given by, referring to (2.26) 


t(y)~ C f ’ (y) ] 1 ^ 2 t c A £ [f (y) ] + dB t ( f (y) ] } 


(2.28) 


Suppose that the condition (2.23) holds. Then, in view of (2.21), the 
eigenvalue problem (2.18) with non-ccnstant U has the asymptotic solution 

* - (l->* 2 )" 1/4 Sinkn/ y [p 2 (s)-l] 1/2 ds, for large n, 
n n n v n 


(2.29) 


where y, * [U(y) - uu/k ] , and k are solutions of the equation 
n n n 


J* [k 2 [l + 2k(-^)U(y) - (^) 2 } 1/2 dy - +nrr i . 


(2.30) 
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We note that (2.30) reduces to (2.19) when U Is constant. For large k , |k| >> u> » 
(2.30) has no real solutions. Thus all higher nodes are improper for a finite 
UU . However, if u> is so large that u> » |k| , (2.30) gives 


or 


J (u) - kU(y)}dy ~ +nTT0r 

-b 


k — {ou 




U(y)dy . 


(2.31) 


which is valid if uu ~ (^^) and n is large. This value of k gives rise to 


P ~*n * 


“ i( ^U L)x + i U L(x * U o t) 
o o 


(2.32) 


1 

where U q • — J U(y)dy is the mean-velocity of the flow. Here we have tacitly 
assumed, by continuation, (2.30) holds also for smaller values of k . The result 
(2.32) shows that at high frequency, the disturbance is a spatially modulated wave 
vith the propagation spped U Q • On the other hand, if |k| , (2.30) has the 

approximation 


k tl • <;> 2 ) 1/2 ^ + (1 . ( u;2)i/2 * - • 


-b 


(2.33) 


which corresponds to the non-propagating modes. 

In the presence of a turning point, we should apply the B.C.'s to (2.28). 

But the resulting characteristic equation is not readily solvable to get 
explicit results as before. 

We remark that, for simplicity, we impose the B.C. p - 0 at y » +b . 

More realistically the field in |y| < b must be matched to that in |y| > b 

at y * + b . This correction is straight-forward but complicates the characteristic 

equation whose explicit solution becomes inaccessible. 

ORIGINAL ^ 
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Instead of simple modes, consider the superposition of modes excited by 
various frequencies 

p(t,x,y) • £ P (e,x,y) , (2.34) 

n 

where 

. i{k («j)x -ujt} 

P n (t,x,y) r j!* a (y,k n .«») t n (y*l« n ^))^ n (aj)« d« (2.35) 

is the n*th mode of p , end k (uu) is the n-th eigenvalue of the BVP (2.12) - (2.13) 

n 

corresponding to the frequency uj . It is easy to verify that ? n given by (2.35) 
satisfies (2.1) after u,v are eliminated. The solution (2.34) may describe the 
disturbance due to a source located at x » 0 . 

Since we are interested in the event near the wave front in the large 
distance down-stream, the asymptotic evaluation of (2.33) will be carried out as 
x -* m , t -» • with (x/t) fixe-*. Let 

t/x - T . (2.36) 

Then (2.33) may be '••ritten as 

-x{k (uu) -urr) 

P n " J* ?_(y. k n .«i)* duu . (2.37) 

n * n n 

where 

• «- 38 > 

Hence, as x -* • , we may apply the method of stationary phase for asymptotic 
evaluation. The stationary points of the phase are determined by the equation 

k' (ui) - t . (2.39) 

ft 

Let ( u* be a solution. Then the contribution to p^ due to this stationary 
point is (15] 
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P„ ~ p„(k*,uu*) [ v i.^,..*TTl 1/2 exp i (k*x - uj*t - 7 Sgnk"(u,*)} 


n n' x|k”Cou* 


where k* ■ k (ui*) . 
n n 

For example, we consider the case of constant U . By (2.19), we have 
k n = " M( a ) ± L n 2 ’ A n = <S )2 " O-M 2 )^) 2 • 


where M - V/cr is the Mach number. It follows that 


V “ 2 1/2 


k" (id) - + a' 2 i' 1/2 (l- -f-) . 
n “ n z 

“ 4 n 


Then the equation (2.39) has the solutions 


id* - +er(arT +M)(1 - M 2 ) (~) [1 + (err +M) 2 ] 1/2 

n ~~ zb 


k * — tu* + 7 ~ . u \ } with T = t/ x 
n n L of - of(on +M) 


Thus we see that, due to boundary absorption, each mode excited by a source 

1 /n 

at x*0 will evolve into space-time modulated waves which decay like |x[ 
downstream. Their respective speeds, referring to (2.42) and (2.43) are given by 


c « |u)*/k | = |M + (M + at/x)’ 1 ! or 


For a variable profile U , if the variation is slow, the WKB method can 
be shown to apply. Thus, if the dispersion equation (2.30) may be solved 
approximately, a similar procedure can be applied to give an asymptotic 


evaluation. 



3. Radiation in an Unbounded planar Jet. 
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In the presence of external source, the FHS of the first equation in (2.1) 
is q(t,x,y) ^ 0 . Let 

f = 1 o XT f (t»x)e~ i ^ kX ~ U)t ^dxdt . (3.1) 

(2tt; 

Then a Fourier transform of (2.1) with q ^ 0 yields (2.4) with the RHS of the 
first equation equal to q t 0 . The eliminated system corresponding to (2.10) 
reads 

L(k)ijr = f - k 2 (l-p 2 H - r* = -ikq , -« < y < “ (3.2) 

R\Jf — 0 at y » + oo (3.3) 


where the B.C. is introduced by the radiation condition for p . Suppose 
G(y,y',k,u0 denotes the Green's function for the BVP (3.2), (3.3). Then 
we have 

f(y,k,tl)) = J’(-ik)G(y,y' ,k,uo)q(y' ,k,tju)dy' . (3.4) 

It follows that 


p(t,x,y) = (ip^’^/J'kTKy.k.uOGty.y 1 ,k,u)) 

X q(y',k,uo)e i(kx ‘ ajt) dy’dkd U ) . 


(3.5) 


By the theory of contour integration, the major contribution to the integral 
(3.5) comes from the singularity set of G . Let be two complementary 

solutions of (3.2) satisfying the B.C.'s at y»os,-« , respectively. Then we 
have 

/^A“ 1 (y')ilr 1 (y)i)r' (y') , y * y' 

G * > , (3.6) 

( k -A _1 (y * ' > » y 5 y' 


12 


ORIGINAL PASS FS 
OF POOR QUALITY 


where A * (t]^2 “ ^2^1^ t * ie Wrons ^^ an w ^ich should be a constant. 

As an example, suppose q is a point source at x = 0, y*0 , so that 
q * q Q (t)6(x)6(y) • Then (3.5) becomes 


P(t*x,y) * (ip^" J7kTKy,k,tjj)G(y,o,k,tu)q o ((u) 

i(kx-tut) 


X e 


dkdu) 


(3.7) 


For a slowly varying velocity U at low Mach number, the result (2.20) by WKB 
method gives 


G(y,o,k,u>) = G(-y,o,k,u>) 


_L_ 2 „ 

2ikg Q g(y) * ’ y ° ’ 


(3.8) 


2 1/2 

where g(y) = [p. (y)-l] , g Q = g(o) . Upon using (3.8) in (3.7), it yields 


p(t,x,y) ~ (-p ) 1 JJ g 1 (k,(jj)Tl(y,k,u)) 


g Q (k,u)) 


g(y,k,uu) 


1/2 


q o (<D> 


ikr y g(s,k,tt))ds +i(kx-uut) 


X e 


dujdk, y i 0 , 


It is difficult to compute the integral in closed form. Sometimes it may be 
possible to evaluate it approximately. 
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4, Propagation in a Slowly Divergent Planar Jet, 

Let the mean-flow be divergent so that, in (1.1), (1.2), 

U * (U,v,0) , 
u ■ (u,v,o) , 


(4.1) 


where U,V are functions of x,y . As a function of x , they are slowly varying. 
At low Mach number, the mean flow is assumed to be incompressible so that there 
exists a stream function j>(|,y) , where § * ex is a slow variable with small 
e > 0 . By definition we have 


u ■ (S.y) B Ms.y) 


y 2 

V - 4 x <5.y> - -4 5 <s.y) - -t^(S,y) • 

In terms of , the governing equations (1.1) - (1.2) can be written in 
components as follows 


(4.2) 


D c u + (4 12 u + ^22 v) 


1 fVp 


p„ d* ’ 


V - (,2 Ki u + , t , i2 v) * - r" ^ ■ 


(4.3) 


TV . /Phi . 3v* 

V + + ■ « 


whe« ■>« ■ <3J + * 2 & - & ' 


I 

l 


When the skew parameter e=0 , the system (4.3) reduces to the parallel flow 
case (2.1). Here we wish to determine the effect of divergence on the 
propagating waves in the jet. By the two-variable method [16], we write 




Then (4.3) may be put in the form 




where 
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D F + B F - Q , 

c c 


0 

dx 

YPod7 

1 _s_ 

p c 5x 

4l2 

^22 

11 
0 o Sy 

-« 2 i 
* m 

-4l2 


q 

o 

o 


First consider the homogeneous B.V.P. in free-space 


D F + B F = 0 , 
6 6 


on which we impose the radiation conditions at “ . 

Surpressing the dependence on t and y , let 

F ■ Z(x,§,e), f * ex , 

s o that 

az = sz az 

dx dx d§ 

Expand Z in a power series in e to get 

Z(x,§,e) - Z o (x,§) + eZ^x.S) + 

In view of (4.8) - (4.10), (4.7) yields the successive approximations 


• . • • 


(4.4) 


<4. 5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 
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(D+B)Z Q - 0 , 

(D+B)Z 1 - Q x ■ -(D 1 +B 1 )Z 0 , 


D - 



» 

D l = 

h-tf 

• 



0 

a 

YP oa^ 

YP o3y 

B = 

-LJL 

p 3x 
K o 

0 

^22 


11. 




_po a y 

U 

u 


0 


0 

B l“ 

i a 
o o a§ 

^*12 

0 


0 

0 

4 12 


(4.11) 

(4.12) 


(4.13) 


Note that (4.11) is the parallel flow problem (2.1). To proceed we consider 
the adjoint problem to (4.11) 


(D* + B*)F* = 0 (4.14) 

where D* and B* are adjoints to D and B , respectively. Now, let F q (x) 
be a solution of (4.11). Then so is 

Z 0 (x,§) «<7(§)F o (x) , (4.15) 

where d7(§) is as yet to be determined. A substitution of (4.15) into (4.12) 
results in 
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(D + B)Z 1 - -(B 1 +D 1 )i7(5)F o (x) . (4.16) 

Let F,G be two vector-functions of t,x,y * Then we define 

<F,G> - J w rf(F'G)dtdxdy , 

where F»G is their dot product. In view of (4.14), we have 

(F*,(D + B)Zj) * <(D* + B*)F*,Z 1 > - 0 . (4.17) 

Thus (4.16) implies that 

<F*,(B 1 + D 1 )F o >C7(§) - 0 (4.18) 

which, noting (4.13), may be written in the form 

a(?) df + S( 5 ) < 7 * 0 ‘ (4 * 19) 

This is the envelope equation describing the slow modulation of waves due to the 
divergence effect. 

Alternatively the problem can be treated by first applying a Fourier transform 
to the systems (4.11) and (4.12). The resulting equations may be reduced to scalar 
equations for the pressure field as done before. This is particularly convenient 
when dealing with the radiation problem (q f* 0) . 

The evelope equation (4.19) is difficult to solve analytically. But it is 
amenable to numerical treatment. For example, in the study of the boundary layer 
stability, Saric and Nayfeh [17] applied two-variable method to account for the 
slow divergent flow. They derived an evelope equation similar to (4.19). Their 
numerical procedures for solution is applicable to the present case. 
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5. Acoustic Waves in an Axisynanetrical Jet. 


ORIGINAL PAGE IS 

OF POOR QUALITY 



The analytical techniques for planar jets may be generalized to an 
axisymmetrical jet. Introduce the cylindrical coordinates (x,r,0) , with the 
corresponding velocity components (U,V,W) and (u,v,w) . 

First consider the cylindrical jet, i.e. U * U(r), V ■ W h 0 . In 
components, the system (1.1) - (1.2) reads 


n j. rdu 1 3 / \ . 1 Sw-i 

Bp + + r S (rv) + t 


Du + — zp + - » T-(rU) ’ 0 , 
p ox r dr 
o 

+ £<■»> - ° • 

Dw + — - 7% - 0 . 

p e r ae 


where D - (£ + u£) . 

For the homogeneous case, q * 0 , we set 


. . i(kx + m0-out) 

(p , u , v ,w) = (p,u,v,w)e 


(5.1) 


(5.2) 


Then (5.1) is reduced to 

ikT]p + yp o {iku + ^ j^(rv) + ^ w} =0 , 
ikTtu + (ik/p Q )p + v(^ j^U) - 0 , 
ikTlv + (l/p o )~ gj^-(rp) * 0 » 
ikTjw + (im/p o r)p * 0 . 

After eliminating u,v,w, in favor of p , the system (5.3) yields 


(5.3) 


i 


* 


* 



4 


<4 
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ORIGINAL PAGL’ iS 

OF POOR QUALITY 


M m (k)? - TlD r Cn" 1 D r )p - a(r)D r p - + k 2 (l-n 2 )]p - 0 


(5.4) 


where 


D r f * r 3r (r ^ ) * ff(r) “ < D r ‘n>/‘H 


(5.5) 


The solutions of (5.4) for m*0,l, known as the axisymmetric and helical modes, 
are of special interest. In contrast with (2.10), Eq. (5.4) is naturally associated 
with an equation of Bessel's type. This is true when U is constant. Then (5.4) 
becomes a Bessel-like equation, but not exactly. 

It is possible to generalize the analytical results in $2 to the present case. 
By (5.5), the equation (5.4) may be written as 


p" + Or " 2T|” 1 ti * )p * - + k 2 (i-n 2 ) +^-}p - o . 

r 1 


(5.6) 


p - Tjijr 


in (5.6) to get 


L m <10* 


(5.7) 


2 

r - {k 2 (l-u 2 ) + + g(r)i); - 0 , (5.8) 


i(r) - - -M . 


(5.9) 


For instance, for large k , we can apply the WKB method or a two-variable 
expansion. The procedure is similar to what was described in §2. One possible 
approach is to seek a solution of the form 


* - a(r)H J [ 1) (r) or 3(r)H< 2) (r) , 


(5.10) 


where \ * yin +1 is the order of the Henkel's functions of first or second 
kind. By substituting (5.10) into (5.8), we may derive equations for a or g , 
to which we apply the asymptotic method. 
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OF POOR QUALITY 

For a divergent jet, the mean flow has the velocity components 

u "J &< r ^» v * W ■ 0 , 


(5.11) 


where |>(g,r), § = ex , is the slowly-varying stream function. The governing 

equations (1.1) - (1.2) for this case takes the form 


where 


and 


V + VP 0 (^ + D r V + i||) -q , 

V + £ + 4 12 u + t 2 2 V " 0 * 

V + ^" V " e2 ^ll u ” *^12 V * 0 

D w + ~ - r§ • 0 , 

« r ae • 

O 

D e M ^ + ^* , 

*1 “If* ^2 V* and 


(5.12) 


(5.13) 


so on. 


Ag.l„ we introduce two variable. « and 5 - « . Thee a two-variable expansion 
procedure a. presented in the previous section may be applied to obtain 
approximate solution to the homogeneous form of (5.12), (q * 0) 


an 
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"Reconstruction of the Mutual Coherence Function for a Moving Source" 
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